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ABSTRACT 


This thesis introduces an integer linear program called the Minimum Cost Flow with 
Congestion Assignment (MCF-CA) model. MCF-CA is a multi-period evacuation model 
that uses a novel approach called congestion assignment to analyze clearing times during 
mass evacuations. Congestion assignment discretizes the nonlinear relationship between 
the number of vehicles on a road segment and the maximum speed at which those 
vehicles can travel. MCF-CA selects among three congestion levels (none, moderate, and 
high) for each road segment in each time epoch. Depending on the congestion level 
selected, MCF-CA limits the number of vehicles that are able to traverse the road 
segment and uses Ak 9 elik’s Time-Dependent Speed-Flow Function (Ak 9 elik 2003) to 
determine the average travel speed of the vehicles for that time period. As a result, we are 
able to determine approximate evacuation clearing times under nonlinear congestion 
effects by solving an integer linear program. We limit residents’ prior knowledge of 
traffic conditions by implementing MCF-CA in a rolling horizon fashion and study the 
impact of this limited knowledge on evacuation patterns. We also model the impact of 
sub-optimal routing decisions on the part of residents by artificially shifting residents 
toward their own shortest paths rather than a “socially optimal” route. 

We find that a mass evacuation can more than double the clearing times of 
individual county evacuations. However, during both county and mass evacuations, 
resident routing choices significantly impact clearing times. As more residents choose 
suboptimal routes, clearing times are prolonged. Lastly, we find that more than 50% of 
residents will experience congestion at some point during the evacuation horizon. 
However, allowing some congestion improves evacuation clearing times by 20-36% over 
not congesting. Although congestion decreases vehicle travel speed by 70-80%, over 
50% more residents are able to start or continue evacuating during each time epoch. 
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EXECUTIVE SUMMARY 


This thesis introduces an integer linear program called the Minimum Cost Flow with 
Congestion Assignment (MCF-CA) model to analyze the impact of mass evacuations on 
clearing times, explore the impact of resident routing choice on clearing times, and assess 
the impact of congestion during a multi-county evacuation. MCF-CA is a multi-period 
evacuation model that uses a novel approach called congestion assignment to analyze 
clearing times. Congestion assignment discretizes the nonlinear relationship between the 
number of vehicles on a road segment and the maximum speed at which those vehicles 
can travel. 

We model a multi-commodity transportation network involving the 216 census 
tracts of Mobile County, Alabama, Baldwin County, Alabama, and Escambia County, 
Florida. We limit residents’ prior knowledge of traffic conditions by implementing MCF- 
CA in a rolling horizon fashion and study the impact of this limited knowledge on 
evacuation patterns. We also model the impact of sub-optimal routing decisions on the 
part of residents by artificially shifting residents toward their own shortest paths rather 
than a “socially optimal” route. 

We find that a mass evacuation can more than double the clearing times of 
individual county evacuations. Exit locations, not simply the number of exits, play a 
significant factor in estimating clearing times. Emergency managers must sometimes 
encourage residents to use exits that are distant in order to ensure timely evacuation of all 
residents. 

We find that routing choices greatly impact clearing times. We develop four 
scenarios that vary the percentage of residents that cooperate towards obtaining a 
“socially optimal” routing or do not cooperate and follow sub-optimal routing. In the base 
case. Global Cooperation, all residents cooperate to obtain optimal routes. In the second 
scenario. Moderate Cooperation, approximately 25% of residents follow their own 
shortest path. In the third scenario. Low Cooperation, approximately 75% of residents 
follow their own shortest path. In the fourth scenario. No Cooperation, 100% of residents 


XV 



follow their own shortest path. Global Cooperation offers the most flexibility in route 
choices. No Cooperation offers no flexibility in route choices. We show that in scenarios 
with a higher percentage of residents following their shortest path, clearing times 
lengthen by as much as two to four times. This result implies that increased flexibility in 
route choice improves clearing time. Additionally, the more flexibility residents have in 
choosing their routes, the less congestion they are likely to experience. Thus, intelligent 
routing decisions by residents are the key to ensuring an efficient evacuation. 

In the best- and worst-case scenarios, more than 50% of residents will experience 
congestion at some point during the evacuation horizon. However, congestion alone does 
not increase clearing times. In fact, if MCF-CA is restricted to disallow congestion, 
clearing time increases significantly. Specifically, allowing some congestion improves 
evacuation clearing times by 20-36%. Although congestion decreases vehicle travel 
speed by 70-80%, over 50% more residents are able to start or continue evacuating 
during each time epoch. 
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I. 


INTRODUCTION 


In 2006, the Federal Highway Administration (FHWA) reviewed and assessed 
63 Federal and State evacuation plans for catastrophic hurricanes and other events 
impacting the states bordering the Gulf of Mexico. The devastating destruction of 
Hurricane Katrina prompted a focused study on the critical issues surrounding mass 
evacuations. The FHWA study, “Catastrophic Hurricane Evacuation Plan Evaluation: A 
Report to Congress,'' assessed the various plans in light of lessons learned about 
evacuation efforts employed during Hurricane Katrina (FHWA 2006). Prior to Hurricane 
Katrina, several mass evacuation plans existed. However, none of them were adequate for 
the scale of destruction that Hurricane Katrina thrust upon Louisiana, Mississippi and 
Alabama (FHWA 2006). Before Katrina made landfall, emergency managers ordered 
more than 1.2 million people to evacuate (DesRoches 2006). 

The FHWA study evaluated the strengths and weaknesses of written evacuation 
plans in terms of how well they met existing Federal Emergency Management Agency 
(FEMA) planning guidelines of the National Response Eramework (formerly National 
Response Plan) (FEMA 2008) and the State and Eocal Guide 101 (EEMA 1996). A 
specific focus of the study was the extent to which evacuation plans were coordinated 
with those of neighboring states and adjoining jurisdictions. The study noted that local 
mutual-aid agreements for joining resources and coordinating decision-making across 
jurisdictions did not address catastrophic events when neighboring jurisdictions would be 
inundated and unable to provide assistance. The study also cited inadequate real-time 
command and control of mass evacuations across all levels of government and across 
multiple states, and it recommended development of additional evacuation modeling 
tools to predict evacuation times and manage operations in real time (EHWA 2006). 

The Catastrophic Hurricane Evacuation Plan Evaluation focused on the written 
contents of evacuation plans, not the quality of the plans in terms of effectiveness and 
efficiency. In 2010, the EHWA published the study “Highway Evacuations in Selected 
Metropolitan Areas: Assessment of Impediment" as a complement to its 2006 study on 

evacuation plans. This study assessed mass evacuation plans for high-threat metropolitan 
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areas to identify and evaluate deficiencies of the national highway system that could 
impede evacuations. The FHWA reviewed existing plans and conducted interviews with 
local jurisdictions, including FHWA Division staff, state and local transportation 
officials, and emergency managers. The expert knowledge of local authorities constituted 
the basis of the report’s findings for each metropolitan area. Results from each area vary 
in the range of impediments described. However, local experts all cited congestion from 
insufficient capacity as the leading impediment to timely evacuation (FHWA 2010). 

The two FHWA studies highlight the impact of road capacity, traffic congestion 
and inter-jurisdiction cooperation in the timely evacuation of large masses of people. 
They also highlight the need for tools to assess congestion and predict evacuation times. 
Effective evacuation modeling tools are essential for state and local emergency managers 
to test scenarios and to facilitate decision-making. Many transportation analysts have 
created evacuation models with the goal of accurately predicting the time required to 
completely evacuate residents, known as the clearing time. We describe four such models 
in Chapter II. However, none of the models focuses on the dynamic complexity of traffic 
congestion at the neighborhood level. 

A. CONTRIBUTIONS 

This thesis introduces an integer linear program called the Minimum Cost Flow 
with Congestion Assignment (MCF-CA) model. MCF-CA is a multi-period evacuation 
model that uses a novel approach called congestion assignment to analyze clearing times 
during mass evacuations. Congestion assignment discretizes the nonlinear relationship 
between the number of vehicles on a road segment and the maximum speed at which 
those vehicles can travel. MCF-CA selects among three congestion levels (none, 
moderate, and high) for each road segment in each time epoch. Depending on the 
congestion level selected, MCF-CA limits the number of vehicles that are able to traverse 
the road segment and uses Ak 9 elik’s Time-Dependent Speed-Flow Function (Ak 9 elik 
2003) to determine the average travel speed of the vehicles for that time period. As a 
result, we are able to determine best-case evacuation clearing times under nonlinear 
congestion effects by solving an integer linear program. We model a multi-commodity 
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transportation system involving the 216 census tracts in Mobile County, Alabama, 
Baldwin County, Alabama, and Escambia County, Florida. We limit residents’ prior 
knowledge of traffic conditions by implementing MCF-CA in a rolling horizon fashion. 
We also model the impact of sub-optimal routing decisions on the part of residents by 
artificially shifting residents toward their own shortest paths rather than a “socially 
optimal” route. 

B. SCOPE OF STUDY 

This thesis conducts a meso-scopic (neighborhood-level) transportation analysis 
of coastal Alabama and Escambia County, Florida. We construct an evacuation road 
network and estimate network vehicle demand using the 2010 Census Bureau American 
Community Survey. We model a road network consisting primarily of freeway segments 
and do not consider interaction of vehicles at segment junctions. We calculate clearing 
times and examine residents’ congestion levels under a variety of scenarios. 

C. THESIS OVERVIEW AND ORGANIZATION 

Chapter II gives an overview of evacuation concepts, reviews existing evacuation 
models, and summarizes basic transportation terminology and traffic flow concepts. 
Chapter III discusses the methodology for developing the road network used in this 
thesis. Chapter IV introduces the MCF-CA model. Chapter V presents detailed results 
and analysis, and Chapter VI presents conclusions and identifies future work. 
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II. BACKGROUND AND LITERATURE REVIEW 


According to a 2004 Cambridge Systematics study, since 1982 both big and small 
cities indicate that they have not been able to keep pace with the rising demand for 
transportation infrastructure (Cambridge Systematics, Inc. 2004). With a growing 
population more dependent on vehicle travel more than at any other time in history, U.S. 
highways, bridges, and roads see more congestion at all times of the day. During normal 
rush-hour traffic, drivers could spend as much as two hours on a 10-mile commute. Even 
during off-peak hours, commuters can spend close to 40 minutes for the same 10-mile 
drive (Cambridge Systematics, Inc. 2004). 

Congestion is a complex phenomenon characterized by stop-and-go or slow-and- 
go traffic. Several conditions directly and indirectly influence congestion intensity, such 
as time of day, lane width, and the presence of car accidents and work zones. In short, 
bottlenecks in the road system create congestion. Bottlenecks form when: 

• Vehicle arrival flow rate exceeds the capacity of a road segment, 

• A queue from a prior bottleneck on the segment has not dissipated, or 

• Traffic flow is affected by downstream conditions (Transportation 
Research Board 2010). 

Regardless of the combination of influencing conditions, congestion results from having 
too many cars in the same place at the same time. 

During mass evacuations, entire regions place simultaneous demands on 
transportation systems. A significant number of vehicles move across a road network in a 
short period of time. State and local emergency managers have emergency evacuation 
plans in place to help alleviate traffic congestion and ensure residents are able to evacuate 
in a timely manner. However, mass evacuations present conditions where overwhelming 
traffic congestion can lead to unsafe road conditions. Vehicle speed decreases rapidly as 
conditions become more congested. Consequently, travel time grows rapidly. Traffic 
congestion can lead to an increase in traffic accidents because of closer vehicle spacing or 
more stalled cars due to overheating during longer commute times (Cambridge 
Systematics, Inc. 2005). 
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This thesis builds upon and complements prior work in network flow 
optimization, traffic congestion analysis, and road capacity concepts. We now review the 
literature most relevant to this study. 

A. EVACUATION MODELS 

Langford (2010) presents a space-time flow optimization model to evaluate 
clearing times for neighborhood-level evacuations. The model determines “best case” 
evacuation routes and clearing times using a minimum cost network replicated over 
several time periods. For the neighborhood analyzed, Langford concludes that “the 
presence of background traffic flow on a major evacuation road with non-evacuation 
traffic does not greatly impact the neighborhood evacuation; rather the overall evacuation 
time is more significantly impacted by the interior roads of the neighborhood” (Langford 
2010 ). 

Yuhas (2011) uses a single-commodity, multi-period minimum cost network flow 
model to evaluate evacuation routes and clearing times for residents of Yolo, Sacramento, 
San Joaquin and Stanislaus Counties of Central California. The model assumes travelers 
have “perfect knowledge” of road conditions before evacuation. Yuhas uses pre-defined 
road capacities to determine the effect of congestion on the network, and he assesses 
contra-flow opportunities to improve clearing times and the effect of highway 
inundations on clearing times. Contra-flow is the reversal of traffic flow on portions of a 
road to allow more flow in the opposite direction. 

Fosgerau (2008) uses a bottleneck model developed by William Vickrey (1969) 
and builds on the bottleneck work of Arnott et al. (1999) who analyzed stochastic 
capacity and demand on congestible roadways. Fosgerau develops explicit expressions 
for the expected marginal and total costs associated with random capacity and demand. 
He concludes that stochastic capacity and demand increases congestion costs up to 50% 
over deterministic capacity and demand (Fosgerau 2008). 

Sheffi et al. (1981) describe the Network Emergency Evacuation Model 
(NETVACl), a macro-traffic simulation model that estimates clearing time during the 
evacuation of the area around a nuclear power plant. The NETVACl simulates a 
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transportation network over several time intervals and uses the mathematical 
relationships between traffic flow-rate (vehicles/hour), speed (miles/hour (mph)), and 
density (vehicles/mile) to evaluate the evacuation process. NETVACl calculates roadway 
capacities according to the Highway Research Board Highway Capacity Manual 
guidelines (Highway Research Board 1965) and updates capacities at the beginning of 
each simulation interval (Sheffi et al. 1981). 

B. ASSESSING CLEARING TIMES 

During hurricanes, state and local emergency managers must ensure complete 
evacuation of vulnerable regions before the onset of pre-landfall hazards such as gale 
force winds (Galveston District, U.S. Army Corps of Engineers 2010). Vulnerable 
regions are those areas that are susceptible to storm surge and flooding. States and 
counties conduct Hurricane Evacuation Studies (HESs) to identify vulnerable regions for 
each hurricane category, to assign each region an evacuation priority, and to estimate 
evacuation clearing times. Accurately assessing the clearing time is essential for safe 
evacuation of residents. 

Emergency managers either mandate or encourage regions to evacuate based on 
the category of the impending storm. The Saffir-Simpson Hurricane Wind Scale 
(SSHWS), originally developed by wind engineer Herb Saffir and meteorologist Bob 
Simpson, lists five categories of storms based on wind intensity. The SSHWS is used to 
assess property damage and the likelihood of injury and death to persons and livestock 
for each category of storm (Saffir-Simpson Team 2012). The time required to clear 
residents from vulnerable regions is a function of the number of evacuating, the 
discretionary decisions of individual residents about when to leave and where to go, and 
the road network’s capacity. 

The decisions of individual residents about when and where to evacuate has been 
studied and modeled by many scholars, but predictable behavior has not been identified. 
Jeffrey Czajkowski developed a multi-period model in which households compare the 
costs of evacuating versus the expected costs of not evacuating in each National 
Hurricane Center forecast advisory period (Czajkowski 2011). Jason Crews also 
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considered a multi-period model and used dynamic programming to develop a generic 
multi-period disaster model in which households decide whether to “stay” or “evacuate” 
in each period (Crews 2012). Dombroski et al. (2006) studied public compliance with 
evacuation orders. The study explored social science factors that could potentially affect 
the behavior and decisions of residents when given the order to evacuate. The general 
consensus of these studies is that more research needs to be conducted to better 
understand the impact of resident evacuation decisions on clearing times. 

The third component impacting clearing time is the road network capacity, which 
is defined as the maximum sustainable hourly flow rate at which vehicles can be expected 
to traverse a road segment in good weather and visibility, with no incidents or accidents, 
no work zone activity, and no pavement deterioration serious enough to affect traffic 
(Transportation Research Board 2010). Roadways operating below capacity are said to be 
under-saturated, while roadways operating above capacity are over-saturated. As 
saturation levels increase, the speed at which vehicles on the road can travel drastically 
decreases. The highly unpredictable routing decisions of individual residents can create 
chaotic traffic patterns which increase saturation levels and, consequently, increase 
clearing times (Fosgerau 2008). 

C. CLEARING TIME STUDIES 

The U.S. Army Corps of Engineers (USACE) conducts traffic flow analyses for 
counties and states to assess the ability of transportation systems to evacuate residents. In 
2010, USACE’s Mobile County District updated its estimated hurricane evacuation 
clearance times for coastal Alabama (Mobile County and Baldwin County) using the 
2010 Census Survey data. Table 1 shows estimated clearing times for five evacuation 
scenarios: one each for Mobile County and Baldwin County evacuating alone; one for 
mass evacuation of Mobile County, Baldwin County and Northwest Elorida; and two 
scenarios for mass evacuation of Mobile County, Baldwin County, Northwest Elorida, 
Mississippi and Eouisiana evacuating with and without T65 contra-flow (Mobile District, 
U.S. Army Corps of Engineers 2010). 
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usage’s evacuation scenarios consider tourist occupancy levels, county 
evacuation zones for various storm categories and the rate at which residents begin to 
evacuate. The road network only includes the primary evacuation route established by 
counties. USAGE researchers estimate and model some non-evacuation background 
traffic but assume local roads can handle most of it. They divide county evacuation zones 
into traffic evacuation zones (TEZs) and estimate the percentage of residents going to 
various destinations. They assign a specific set of routes to each TEZ based on the 
destination of its residents. Glearing times ranged from 10 to 47 hours. The study 
determined that the bottleneck for evacuating both Mobile Gounty and Baldwin Gounty is 
along Interstate Highway 65 at Highway 113 in the neighboring Escambia Gounty, 
Alabama. Table 1 identifies other critical road segments during each scenario. 


Clearance Times Summary 



Critical 

Roadway Segment 

Call 

Range 

Cat 2 

Range 

Cat 3 

Range 

Cat 4 

Range 

Cats 

Range 

Baldwin 

County 

Ala 59 thru Loxley 

Includes Mobile thru-traffic 

No Interstate impacts 

10 to 22 

hours 

15 to 29 

hours 

15 to 29 

hours 

22 to 36 

hours 

22 to 36 

hours 

Mobile 

County 

1-10 eastbound to Baldwin Co. 

Includes Baldwin thru-traffic 

No Interstate impacts 

15 to 20 

hours 

15 to 20 

hours 

17 to 23 

hours 

20 to 25 
hours 

20 to 25 

hours 

Additional 

States 

1-65 eastbound with FL w 1-10 eastbound 

(worst case) 

Includes LA, MS (eastbound) and FL (westbound) 

15 to 25 

hours 

21 to 32 

hours 

25 to 36 

hours 

34 to 45 

hours 

36 to 47 

hours 

1-65 eastbound with FL 

(w/o 1-10 eastbound) 

Includes only FL (westbound), no LA or MS (eastbound) 

15 to 25 

hours 

21 to 31 

hours 

24 to 35 

hours 

33 to 44 

hours 

34 to 45 

hours 

Contraflow 

1-65 eastbound with FL w 1-10 eastbound 

Includes LA, MS (eastbound) and FL (westbound) 

15 to 22 

hours 

15 to 29 

hours 

17 to 29 

hours 

22 to 36 

hours 

23 to 36 

hours 


Note: Contraflow benefits on 1-65 will not generate clearance times lower than the highest time at the county critical roadway segment for 
Baldwin or Mobile County. 


Table 1. U.S. Army Gorps of Engineers 2010 clearing time estimates for coastal 
Alabama. Erom Mobile District, U.S. Army Gorps of Engineers, 2010. 


D. ROADWAY AND CAPACITY CONCEPTS 

The Transportation Research Board Highway Capacity Manual (HCM) assigns 
free-flow speed, the speed at which vehicles are unimpeded by other vehicles on the road, 
and capacity to various types of facilities (roadways) (Transportation Research Board 
2010). It also divides them into two general facility types: uninterrupted and interrupted. 
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The frequency of signalized intersections characterizes a facility’s type. Uninterrupted 
facilities are freeways, multilane highways, and two-way highways. Freeways do not 
have signalized intersections, and multilane highways and two-way highways have at 
least two miles between signals. Interrupted facilities consist of urban streets, where 
signalized intersections are less than two miles apart (Transportation Research Board 
2010 ). 

This thesis focuses on uninterrupted facilities. However, approximately 6% of the 
road segments in the counties studied are interrupted facilities. Modeling interrupted 
segments involves methodologies for incorporating various signalized delays, mid-block 
free-flow speeds, lane width adjustments, and other adjustment factors that are beyond 
the scope of this thesis. For simplicity, we model these segments as uninterrupted 
facilities. 


E. TRAFFIC FLOW CONCEPTS 

Traffic flow is a complex phenomenon encompassing the interactions of vehicles, 
bicycles, and pedestrians sharing transportation infrastructure. The key to understanding 
traffic flow and congestion is to understand the relationship between traffic flow rate, 
speed, and density (Hall 1992). Flow rate {q) measures the number of vehicles passing a 
reference point over a period of time. Speed (v) measures the distance traveled over a 
period of time. Density (D) measures the number of vehicles occupying a length of 
roadway in a particular instant (Transportation Research Board 2010). 

Transportation analysts use a variety of traffic flow equations to represent the 
relationship between the flow rate, speed and density. Wardrop (1952) developed the 
fundamental flow-speed-density equation: 

D = ^ (2.1) 

V 


where 
D = 

Q = 


V 


density [vehicles/mile/lane] 
flow rate [vehicles/hour/lane] 
mean speed of traffic [mph]. 
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The HCM utilizes this widely-accepted equation as the basis for its speed-flow 
models. Speed-flow models predict speed as a function of flow rate. They are centered on 
the concept of free-flow speed (Transportation Research Board 2010). The underlying 
assumption of Wardrop’s equation is that vehicles move together harmoniously at 
constant speed. While this is approximately true for under-saturated conditions, over¬ 
saturated conditions violate the equation’s implied assumption that vehicle spacing 
remains constant in over-saturated conditions (Hall 1992). Vehicle spacing in over¬ 
saturated, slow-and-go congested traffic is not constant. Therefore, HCM methodologies 
are inadequate to model over-saturated conditions. 

F. AK^ELIK’S TIME-DEPENDENT SPEED-FLOW FUNCTION 

Rahmi Ak 9 elik developed the Time-Dependent Speed-Flow Function (TDSF), 
which is both a steady-state queuing delay function for under-saturated conditions and a 
deterministic delay function for over-saturated conditions (Ak 9 elik 2003). The TDSF 
estimates travel speed at a given saturation level (called the degree of saturation, x). The 
main parameters describing the TDSF are: 

V = speed at a given degree of saturation v [mph] 

Vf = free-flow speed [mph] 

= speed at capacity [mph] 

= duration of the analysis period [hours] 

Q = nominal capacity [vehicles/hour] 

= degree of saturation below which the traffic delay is zero [dimensionless] 

= delay parameter [vehicles/mile^] 

V, = initial queued demand [vehicles] 

= flow-rate [vehicles/hour] 

V = degree of saturation [dimensionless] 

= (da / 0 ) 

x' = degree of saturation adjusted to account for V,. effects [dimensionless] 

= x + [N,/(QT^)] 

z = a parameter defined for convenience [dimensionless] 
z = v-l + 2V,/(er^) 


(Ak 9 elik 2003). 
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The TDSF is: 


v=Vfl{\ + 0.15VfTf[z + {z^ + %k,{x-xjl {QTf) + \6k,NJ {QT^ff]} 



forx'>x^ 

(2.2) 


forx'Kx^ 

(2.3) 


(Ak 9 elik 2003). 

When there is no initial queued demand (N^ = 0), Equation (2.2) simplifies to: 

V =v^/{I + 0.25 v^7}[z + (z^ + 8k^(x-xj/(e7}))^]} for x'> x^ (2.4) 
(Ak 9 elik 2003). 

Ak 9 elik calibrates this model by using empirical data to assign the appropriate 
values of v„ / and x„. He proposes using the following values: 

“(i) use the same value of speed ratio for all four classes of basic freeway 
segments and for all four classes of multilane highways (v„ / = .85 ) for 

freeways and 0.82 for multilane highways have been selected); and 

(ii) use the same value of degree of saturation to determine the flow limit 
for free-flow speed (or zero traffic delay) for all four classes of basic 
freeway segments and for all four classes of multilane highways (A,= 

0.70 for freeways and 0.65 for multilane highways have been selected).” 

Ak 9 elik also suggests using x^= 0.50 for urban streets. Singh (1995) proposes using the 

values in Table 2 for the delay parameter, : 


Roadway Type 

Free-flow 
Speed (mph) 

kd 

Freeway 

75 

0.1 

Arterial 

(unintermpted) 

62 

0.2 

Arterial 

(intermpted) 

50 

0.4 

Secondary 

(intermpted) 

37 

0.8 

Secondary 
(hi^ friction) 

25 

1.6 


Table 2. Representative Parameters for TDSF. After Singh, 1995. 
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III. MASS EVACUATION NETWORK MODEL 


This chapter presents the development of the road network studied in this thesis. 
We model evacuating residents of the 216 census tracts (neighborhoods) of Mobile 
County, Baldwin County, and Escambia County as outlined in the 2010 Census Survey. 
The network modeled includes interstate highways, major county roads, and minor 
county roads. We assume that local roads can handle non-evacuating background traffic. 

A. BUILDING THE NETWORK 

We start with a collection of road segments extracted from the 2012 U.S. Road 
Data Base for Network Modeling (Brown and Halwachs 2012). The data for each road 
segment includes the segment’s length in miles and its speed limit in miles per hour 
(mph). We use nominal capacities as suggested by HCM (2010) and Singh (1995) 
(detailed in section B of this chapter). Because transportation networks change over time, 
we use Google Earth, a free online geographical information software (Google, Inc. 
2012), to visually verify the existence and number of lanes for each road segment. We 
keep all road segments verified to exist and discard all others. We manually add road 
segments where necessary to accurately represent the current road network. 

We model the road network using a set of nodes and a set of arcs. Source nodes 
have supplies that can be routed through the network to satisfy demand at sink nodes. 
Nodes that have neither supply nor demand are called transshipment nodes. We identify a 
source node for each census tract (see Eigure 1). All residents within a census tract 
originate from the assigned source node. Using each county’s primary evacuation route, 
we designate exit points outside the vulnerable area as sink nodes. We assign to each 
segment of the primary evacuation route the number of lanes visually verified via Google 
Earth. We assign one lane to each segment that is not part of the primary evacuation 
route. 
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Figure 1. Partial Baldwin County Google Earth spatial representation (a) 
and Census Tract map (b). 
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B. ROAD SEGMENT PARAMETERS 


Each road segment has associated parameters that describe its capacity and speed. 
For each road segment, we determine the values of these parameters in under-saturated 
(free-flow) conditions. The HCM recommends free-flow speeds and capacities for 
freeway facilities as shown in Table 3. 


Free-Flow 
Speed (mph) 

Capacity per lane 
(ve hide s/hour) 

70 

2400 

65 

2350 

60 

2300 

50 

2100 


Table 3. HCM freeway facility speeds and capacities. 

After Transportation Research Board, 2010. 


Singh (1995) also recommends free-flow speed and capacities for various facility 
types as shown in Table 4. 


Type 

Free-Flow 
Speed (mph) 

Capacity 
(vehicles/hour) 

Freeway 

75 

2000 

Arterial 

(uninterrupted) 

62 

1800 

Arterial 

(interrupted) 

50 

1200 

Secondary 

(interrupted) 

37 

900 

Secondary 

(high-friction) 

25 

600 


Table 4. Free-flow speeds and capacities for various 
facility types. After Singh, 1995. 


15 




Based on these recommendations, we divide road segments into four categories 
(I, II, III, and IV) based on their nominal speed and assign to each class a nominal 
capacity as shown in Table 5. 


Road Segment 
Category 

Nominal Speed 
(mi*) 

Nominal Capacity per 
lane (vehicles/hour) 

I 

60 + 

2300 

II 

50-59 

1700 

III 

40-49 

900 

IV 

30-39 

600 


Table 5. Road segment classes used in this study. 


We model three congestion levels: no congestion, moderate congestion, and high 
congestion. When a road segment is operating below its saturation level, no congestion 
occurs and there is no traffic delay. Moderate congestion occurs when a segment operates 
at 100% saturation. High congestion occurs in over-saturated conditions; based on 
Ak 9 elik’s TDSF, we use 106% of each segment’s free-flow capacity to define high- 
congestion conditions. Table 6 shows that as congestion levels increase, arc capacity 
increases and the speed at which vehicles travel decreases. 



1 CONGESTION LEVELS | 

1 NONE 1 

1 MODERATE | 

1 HIGH 1 

Type 

Nominal 

Speed 

(mph) 

Capacity 

(vehicles/hour) 

Trawl Speed 
(mi*) 

Capacity 

(whicles/hour) 

Trawl Speed 
(mph) 

Capacity 

(whicles/hour) 

Trawl Speed 
(mph) 

1(a) 

65 

1610 

65 

2300 

28 

2438 

18 

1(b) 

60 

1610 

60 

2300 

27 

2438 

17 

n(a) 

55 

1105 

55 

1700 

19 

1802 

14 

n(b) 

50 

1105 

50 

1700 

18 

1802 

13 

m(a) 

45 

450 

45 

900 

10 

954 

8 

m(b) 

40 

450 

40 

900 

9 

954 

8 

IV 

30 

300 

30 

600 

6 

636 

5 


Table 6. Sample of segment capacity and travel speed for the three congestion levels. 


The moderate congestion level decreases travel speed by more than half of the 
nominal speed. Figure 2 depicts the speed-flow curve using the Ak 9 elik TDSF. It shows 
that travel speed decreases rapidly as saturation levels increase. 
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Demand flow rate, q (veh/h) 


Figure 2. TDSF speed versus flow-rate for basic freeway segments. From Ak§elik, 2003. 


We calculate the time (t) required to traverse a segment as its length (/) divided 
by the travel speed at the given saturation level (v) : 


t = l/v 


( 3 . 1 ) 


C. NETWORK ASSUMPTIONS 

For simplicity, we assume roads are empty of traffic when evacuation starts. 
However, we can model initial traffic by converting the appropriate transshipment nodes 
to source nodes. Additionally, we do not consider highway on-ramp capacity restrictions. 
We can model on-ramp capacity restrictions by using additional arcs with lower capacity 
at the appropriate nodes in the network. We assume that residents evacuate immediately, 
and we model all vehicles as standard passenger vehicles. 
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D. NEIGHBORHOOD SUPPLIES 


We calculate the number of vehicles evacuating from each census tract using the 
2010 Census Survey vehicle availability statistics. The Census Survey divides households 
into four categories: 

• Households with 1 available vehicle, 

• Households with 2 available vehicles, 

• Households with 3 available vehicles, and 

• Households with 4 or more available vehicles. 

For simplification, we assume that all households with 4 or more vehicles have only four 
vehicles available. 

For each household /, the Census Survey provides data on the number of available 
vehicles a^and the number of potential drivers o-. Using this data, we calculate , the 
number of drivable vehicles in household i, as 

di=rmn(a.,Oi) (3.2) 

Let be the number of households in a given census tract with exactly d 
drivable vehicles, where <i = l,2,3,4. We assume that each household with available 
vehicles will use at least one vehicle to evacuate; thus, the minimum possible supply 
originating in the tract is calculated as 

5 . =yH, (3.3) 

mm d ^ ' 

d 

The maximum supply introduced in a tract is simply the total number of drivable 
vehicles in the tract: 

(3.4) 

d 

Households with multiple occupants may choose to evacuate using one vehicle to ensure 
that the family evacuates together, or they may use more than one to afford themselves 
additional flexibility and to remove additional belongings from the vulnerable region. For 
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simplicity, we calculate each neighborhood’s supply V as the average of the 
neighborhood’s minimum and maximum possible supply: 




(3.5) 


Figure 3 provides an example of network demand calculation for a census tract 
with 1420 households with vehicles. 


=1420 

d 


= (1420 + 2730)/2 
= 20"5 


Figure 3. Example supply calculation for a census tract. 


# of 

# Households 

Dnvable 

with exactly d 

Cars 

dnvable cars 

(d) 


1 

458 

2 

688 

3 

200 

4 

74 


E. DISCUSSION 

The evacuation model is a large, multi-commodity road network. The road 
network consists of 2710 segments and 825 intersections (i.e., nodes). Commodities are 
called resident types and represent the 216 census tracts of Mobile County, Baldwin 
County, and Escambia County. In total, the counties evacuate 433,334 residents: 
197,483 for Mobile County, 95,304 for Baldwin County, and 140,547 for Escambia 
County. 

In the next chapter, we develop an evacuation modeling tool that predicts clearing 
times and minimize the total travel time of evacuating residents. 
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IV. MINIMUM COST FLOW WITH CONGESTION 
SCHEDULING (MCE-CA) 


This chapter details the standard Minimum Cost Flow (MCF) model and presents 
the Minimum Cost Flow with Congestion Scheduling (MCF-CA) formulation, which 
models congestion during mass evacuations. 

A. MINIMUM COST FLOW 

The MCF-CA model builds on the standard MCF model and its notation. Thus, 
we briefly review the MCF formulation. The goal of MCF is to minimize the total cost to 
move supplies from source nodes to sink nodes (Ahuja et al. 1993, 4-6) 

I. Notation 

Following the conventions in Ahuja et al. (1993) let G = (N,A) denote a graph, 
where N is the set of nodes, indexed by n (alias i and j), and A is the set of directed arcs 
Each node n is a junction point that connects two or more arcs. An arc is a segment 
that provides a path for flow from one node to another. Directed arcs allow flow in only 
one direction. When flow is allowed in two directions, we use two directed arcs. For 
example, if a road network has traffic flow in eastward (i -»j) and westward (/<-j) 
directions, then we include arc (ij) and arc (j,i) are in the set of directed arcs. 

Let denote the capacity, or maximum allowable flow, on arc (iJ), where 
u.j > 0. Let c,^ denote per-unit cost of flow along arc (iJ). Let denote the flow along arc 
(iJ). Let s. denote the exogenous supply at node i (see Ligure 2). 



Ligure 4. Dumbbell graph representation illustrating flow from node i to node /. 

After Yuhas, 2012. 
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2 . 


Formulation 


The goal of MCF is to minimize cost: 

MCF: 


Z = min V c,y, 

Y y y 


(4.1) 

Subject to 



Zr, - Zr.=^. 

j^N i^N 

N 

(4.2) 

0<Y,<u^ 

V(j,7)eA 

(4.3) 


Constraint (4.2) is a flow balance constraint. For each node n, the total out-flows 
minus in-flows must equal the exogenous supply at node n. Constraint (4.3) bounds the 
flow along each arc (ij) using the arc’s capacity. 

B. MINIMUM COST FLOW WITH CONGESTION SCHEDULING (MCF- 

CA) 

In contrast to MCF, MCF-CA is a multi-period model. While MCF’s objective is 
to minimize the cost to move flow through the network, MCF-CA has two objectives. 
The primary objective is to minimize the number of residents who do not successfully 
evacuate during the planning horizon, and its secondary objective is to minimize the total 
cost to move flow through the network. In an evacuation scenario, time is critical. Thus, 
we use the time required to travel each arc (i,j) as its cost. We model the residents of the 
216 census tracts of Mobile County, Baldwin County and Escambia County as individual 
commodities, and we ensure that for each commodity, the total time spent traveling is 
appropriate given the length of each time epoch. 

In order to model nonlinear congestion effects, MCF-CA selects a congestion 
level for each arc in each time epoch. This congestion level acts to limit the arc’s capacity 
and uses Ak 9 elik’s TDSF Equation (2.4) to determine the average travel speed. MCE-CA 
uses Equation (3.1) to calculate travel time for each arc based on the congestion level 
selected. Eor simplicity, we introduce an artificial node, ‘'safe that is connected to each 
sink node via an arc with an appropriately large capacity to allow all flow to traverse it 

(see Eigure 5). The cost per unit of flow to traverse an artificial sink arc is zero. 
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Figure 5. Network with an artifieial “safe” node. 

1. Formulation 
MCF-CA: 

Index Use [cardinality]: 


riG N 

nodes, alias {i, j) 

[825] 

tgR 

resident types (commodities) 

[216] 

t^T 

time epochs 

[~ 6-130] 

Ig L 

congestion level 

[3] 


directed arc from node i to node j 

[2720] 

Data [units]: 

^mt 

supply of resident type r at node n during epoch t 

[vehicles] 

^ijl 

capacity of arc (/, j) at congestion level 1 

[vehicles / hr] 

^ijl 

time to travel arc (i, j) at congestion level 1 

[minutes] 

Pn 

penalty for being stranded at node n 

[minutes] 

e 

time length of each epoch 

[minutes] 
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Decision Variables [units]: 

flow of resident type r on arc (/, j) at congestion level I [vehicles] 

during epoch t 

number resident type r stranded at node n after epoch t [vehicles] 

D^jij binary (1 if congestion level I is chosen, on arc (/, j) [binary] 

during epoch t; 0 otherwise) 


Formula tion 

Z = min 

Y,B,D 


z 


r Y 

ijl ^ ijrlt 


+ e Y, {^ + Pn)K 


IeLJgT 


nEN,rER, 

IeT 


(4.4) 


Subject to 


II 

T 

+ 5 - Vy + yy 

nrt njrlt inrlt 

{n,j)EA, (/,«)gA, 

IeL leL 

\/n G N,r e R,t gT 

(4.5) 

•IS' 

s 

VI 

^ijit 

\/(i, j) G A,l gLpgT 

(4.6) 

X ^ijlt ^ 

leL 


\/(i,j)GA,tGT 

(4.7) 

ijrlt 

(iJ)^A, 

IeL 

neN 

\frGR,tGT 

(4.8) 

^,.^{0,1} 


\/(i, j) G A,l gLjgT 

(4.9) 

IV 

o 


V(/, y) G A,r G R,l G L,t gT 

(4.10) 

B >0 

nrt 


\/nG N,r G R,t g T 

(4.11) 


2. Objective Function 

MCF-CA’s objective function (4.4) contains two terms. The first term mi rrors the 
objective function of MCF and minimizes residents’ travel time through network. This 
term serves two purposes in MCF-CA. First, it ensures that residents take reasonable 
paths during their evacuation. Second, it helps to ensure that the binary variables are 

set correctly. We elaborate on this function in section IV.B.4.b. 

The second term in (4.4) reflects MCF-CA’s primary objective, which is to ensure 
that residents make as much progress as possible toward exiting the vulnerable area 
during the planning horizon. We assign to each node n a stranding penalty 
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proportional to the cost (in minutes) of the shortest path from n to ‘safe’ in maximum 
congestion. This penalty encourages residents to move closer to a sink node during 
epochs when the network does not have enough capacity for complete evacuation. 
Successfully evacuated residents are stranded at the ‘safe ’ node, with a penalty cost of -2 
to reward residents for successful evacuation. The stranding penalty term is multiplied by 
the epoch length, e, in order to ensure that the primary objective to evacuate all residents 
supersedes the secondary objective of minimizing travel time. The epoch length e is a 
sufficiently large multiplier because residents cannot travel longer than the epoch time 
length in any given period. 

3. Constraints 

Constraint (4.5) is similar to the balance of flow Equation (4.2). It is a 
bookkeeping constraint that records the number of residents stranded at node n at the end 
of each epoch. Constraint (4.6) is similar to constraint (4.3) in that it bounds the flow on 
each arc to the arc’s capacity. Constraint (4.7) selects exactly one congestion level for 
each arc. Constraint (4.8) limits the average cumulative travel time of each resident type 
to at most the epoch length in order to ensure that residents do not travel excessively 
long. This is a relaxed constraint for limiting each vehicle’s travel time. In principle, 
residents could be modeled as 433,334 individual commodities, and each vehicle’s travel 
time could be limited. However, microscopic evacuation patterns are beyond the scope of 
this thesis. Constraints (4.9)-(4.11) declare variables types. 

4. Discussion 

Here we discuss the role of the objective function value and provide an example 
of how the appropriate congestion level is selected in order to minimize clearing time. 
We also discuss model implementation and solution times. 

a. Objective Function Value and Clearing Time Estimation 

MCF-CA’s objective function value does not explicitly calculate clearing 
time. Rather, it is designed to reflect the behavior of residents with imperfect knowledge 
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of future conditions but a desire to move quickly toward safety. It reflects the cumulative 
travel time of all residents and the penalties incurred for stranding residents in each time 
period. 

Given an optimal solution to MCF-CA, we calculate clearing time during 
post-optimality analysis. The approximate clearing time (rounded to the next hour) is 
simply the number of 60-minute epochs required to evacuate all residents. We can derive 
tighter bounds for the clearing time by analyzing the flows in the last epoch; however, 
such analysis is beyond the scope of this thesis. 

b. Selection of Congestion Level 

MCF-CA limits the flow on each arc such that the sum of all flows on the 
arc does not exceed the arc’s capacity, given its congestion level. Note that by itself, 
constraint (4.6) allows MCF-CA to select a higher level of congestion than is required to 
accommodate the optimal flows. For example, let congestion levels 1=1, 1=2, and 1=3 
denote the states of no congestion, moderate congestion, and high congestion, 
respectively. Let F.j, denote the total flow on arc (i, j) in epoch t, and suppose that 

arc ((, j)experiences total flow = 'Yu^ijrit optimal solution. MCF-CA can 

r,l 

satisfy constraint (4.4) by setting either Dy^, = 1 or Dy^, = 1. However, if Fy^ > 0 then 
we must have Dy^,in an optimal solution, since Dy^^ = 1 would result in higher (worse) 

value for the first term in (4.4) and no change in the second term in (4.4). Furthermore, 
the left hand side of constraint (4.8) is greater when Dy^, = 1 than when Dy^, = 1. Thus, 

constraint (4.6) and the first term of the objective function (4.4) work together to ensure 
that the variables Dy,^ accurately reflect the minimum level of congestion required to 

evacuate all residents in minimum time for all for which Fy, > 0 . Since MCF-CA 

can set Dy^^ arbitrarily when Fy^ = 0, we analyze congestion level selections for only 
those combinations for which F . > 0 . 
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c. 


Model Implementation 


We solve MCF-CA with a 1-epoch lookahead using GAMS/CPLEX 12 
and run it on computers equipped with an Intel 3.0GHz processor and 96GB of RAM. 
The largest single-county evacuation, Mobile County, generates over 98,000 constraints, 
470,000 variables, 2,000,000 non-zeroes, and over 3,000 discrete variables for each 60- 
minute time period; its solution time is approximately 12 hours. The mass evacuation 
generates over 189,000 constraints, 1,900,000 variables, 9,000,000 non-zeroes, and over 
8,000 discrete variables for each 60-minute time period; its solution time is 
approximately 2.5 days. 
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V. RESULTS AND ANALYSIS 


We now use MCF-CA to analyze clearing times for Mobile County, Baldwin 
County, and Escambia County during both single-county evacuations and during a three- 
county mass evacuation. 

A. COMPARISON OF COUNTY AND MASS EVACUATIONS 

USACE measures the number and destination of evacuees by comparing normal 
everyday traffic patterns to the traffic patterns during the evacuation, both within the 
region and immediately outside of the region near exit points (Mobile District, U.S. Army 
Corps of Engineers 2010). They attribute traffic volume spikes to the movement of 
residents from the evacuating region. During county-only evacuations, residents compete 
with normal traffic volume just outside the county at exit points. However, during mass 
evacuations of multiple counties, the road segments immediately outside of a county are 
already experiencing increased traffic volume due to the evacuation efforts of the 
neighboring county. This forces residents to compete for access to road segments with 
traffic volumes well above the normal pattern. 

Eor each county, we solve MCE-CA to establish a best-case clearing time when 
residents follow optimal routes and no road impediments occur. Eigure 6 displays the 
clearing times for all counties evacuating in isolation and during a mass evacuation. As 
the figure indicates, Escambia County’s best-case clearing time is 7 hours during a 
single-county evacuation and 11 hours during a mass-evacuation; Mobile County’s best- 
case clearing time is 14 during a single-county evacuation and 15 hours during a mass 
evacuation; and Baldwin County’s best-case clearing time is 6 during a single-county 
evacuation and 15 hours during a mass evacuation. As expected, clearing times for the 
mass evacuation are substantially greater than those of county evacuations. Clearing 
times increase by 7-133% when counties evacuate together. Baldwin County bears the 
brunt of mass evacuation because it is geographically sandwiched between Mobile 
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County and Escambia County. Residents from both Mobile and Escambia Counties flow 
through Baldwin to evacuate. The influx of residents significantly impede the evacuation 
flow of Baldwin County residents. 


Best-Case Clearing Times 
Single-County Evacuation vs. Mass Evacuation 



■ Mass Evacuation 
□ Single-County Evacuation 


Eigure 6. Best-case single-county and mass evacuation clearing times. 


Eigure 7 shows the hourly percentage of residents evacuated during county and 
mass evacuations. By hour 7 during single-county evacuations, both Baldwin County and 
Escambia County are completely evacuated. However, during mass evacuations, these 
counties have evacuated only 70% and 80% of their residents, respectively. By hour 7, 
Mobile County has evacuated only slightly fewer residents during the mass evacuation, 
41% versus 51%. 
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Single-County vs. Mass Evacuation Progression 



Hours Elapsed 


SINGLE-COUNTY EVACUATION 

Escambia 

Baldwin 

Mobile 

MASS EVACUATION 
“““ Escambia 
Baldwin 
“““ Mobile 

All Counties 


Figure 7. Single-county clearing times compared to mass evacuation clearing times. 


Escambia County’s clearing time increases by 57% in part because three of the 
nine exit points utilized during its single-county evacuation flow into Baldwin County 
(see Figure 8). Figure 9 reveals that 39.5% of Escambia County residents use the exits at 
Hwy 98W, Hwy 90W, and I-lOW during single-county evacuations. However, these exits 
are still within the vulnerable region during mass evacuations involving Baldwin County. 
Consequently, residents either reroute to other exits within the county or reroute to the 
exits of the adjacent counties. Residents increase their use of exits within the county by as 
much as 80%. Additionally, approximately 7% use Baldwin County’s I-65N exit. Both 
adjustments result in increased flow through fewer exits, thereby, prolonging the time 
required to move all residents out of the vulnerable area and increasing clearing times. 


31 

















Figure 8. Escambia County’s exits I-lOW, U.S.-90W, and U.S.-98W are located on 

the Baldwin County border and thus do not exit the vulnerable region during 
a mass evacuation. After Google, Inc., 2013. 


Escambia County Exit Use 


County Evacuation 


Mass Evacuation 



■ Escambia: I-lOE 

■ Escambia: I-low 

■ Escambia: US-90E 

■ Escambia: US-90W 

■ Escambia: US-98E 

■ Escambia: US-98W 

■ Escambia: FL-4E 

■ Escambia: US-29N 

■ Escambia: FL-97N 



■ Escambia: I-lOE 

■ Escambia: I-lOW 

■ Escambia: US-90E 

■ Escambia: US-90W 

■ Escambia: US-98E 

■ Escambia: US-98W 

■ Escambia: EL-4E 

■ Escambia: US-29N 

■ Escambia: FL-97N 

■ Baldwin: I-65N 

■ Mobile: US-98W 

■ Mobile: US-45N 

■ Mobile: US-43N 


Escambia Exits 


Escambia Exits 


Figure 9. Escambia County evacuation exit use for all scenarios. 
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Mobile County’s clearing time is minimally impacted during the mass evacuation. 
Although residents lose the original Mobile County I-65N exit just outside the county, 
they take advantage of Baldwin County’s north-bound contra-flow and move 29% of 
residents through Baldwin County I-65N exit (see Figure 10). Mobile County residents 
travel an additional 22 miles to the Baldwin County I-65N exit. Flowever, this shift 
benefits Mobile County residents by alleviating the need to reroute as many residents to 
the other exits within the county which would prolong clearing times, as seen in the case 
of Escambia County. Figure 11 shows the proximity of the two exits. Flere, the close 
proximity of the exits benefits the mass evacuation effort of Mobile County. 

Baldwin County’s clearing time more than doubles during mass evacuations 
because it loses one of its two exits (I-lOW). As a result, an additional 10% of Baldwin 
County residents use the remaining I-65N exit (see Figure 12). The remaining 45% who 
change their exit travel as much as 60 miles through either Mobile County or Escambia 
County to evacuate. 


Mobile County Exit Use 


County Evacuation 



■ Mobile: I-65N 

■ Mobile: US-98W 

■ Mobile: US-45N 

■ Mobile: US-43N 


Mass Evacuation 



■ Mobile: I-65N 

■ Mobile: US-98W 

■ Mobile: US-45N 

■ Mobile: US-43N 

■ Baldwin: I-65N 

■ Escambia: I-lOE 

■ Escambia: US-90E 

■ Escambia: US-98E 

■ Escambia: FL-4E 

■ Escambia: US-29N 

■ Escambia: FL-97N 


Mobile Exits 


Mobiie Exits 


Figure 10. Mobile County exit use during county and mass evacuations. 
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Figure 11. Baldwin County I-65N exit location and Mobile County I-65N exit 
location on the border of Baldwin County. After Google Earth, 2013. 


These findings demonstrate the significance of exit locations during evacuation. 
During mass evacuations, some of the exits that residents rely upon for county 
evacuations are still within the vulnerable region. Consequently, either residents reroute 
to the remaining exits in the county or travel upwards of 60 miles to evacuate through the 
neighboring county. In the case of Mobile County, residents may travel two counties over 
to exit the vulnerable region. In either circumstance, loss of exits during mass 
evacuations increases clearing times substantially. 

We note that this solution is only one of multiple optimal solutions. The solver 
does not discriminate between the various resident types. Therefore, when more than one 
type of resident is stranded at a node in the same epoch, the solver arbitrarily chooses 
which residents utilize which routes. Moreover, the solver could generate a different 
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optimal solution by choosing a different combination of resident types to move. 
Emergency managers must keep in mind that multiple optimal solutions may exist for the 
evacuation problem. 


Baldwin County Exit Use 


County Evacuation 



Baldwin Exits 



■ Baldwin: now 

■ Baldwin: I-65N 

■ Escambia: I-lOE 

■ Escambia: US-90E 

■ Escambia: US-98E 

■ Escambia: FL-4E 

■ Escambia: US-29N 

■ Escambia: FL-97N 

■ Mobile: US-98W 

■ Mobile: US-45N 

■ Mobile: US-43N 


Baldwin Exits 


Figure 12. Baldwin County exit use during county evacuations. 


B. VALIDATION AND COMPARISON WITH EXISTING MODELS 

We compare MCF-CA’s clearing time results for Mobile County and Baldwin 
County single-county evacuations to those described in the USACE 2010 study of coastal 
Alabama (Mobile District, U.S. Army Corps of Engineers 2010). We find that in general, 
MCF-CA’s clearing times tend to be lower than USACE’s clearing times. 

MCF-CA’s best-case clearing time for Mobile County is 14 hours and uses a 
network supply of 197,483 vehicles evacuating to outside of the county. The comparable 
USACE Category 5 Immediate Response scenario generates a 20-hour clearing time 
using a network supply of 132,932 vehicles evacuating to outside of the county (Mobile 
District, U.S. Army Corps of Engineers 2010). The MCF-CA clearing time for Mobile 
County is 30% shorter than USACE’s clearing time. 

Similarly, MCF-CA’s best-case clearing time for Baldwin County is 6 hours and 

uses a network supply of 95,304 vehicles evacuating to outside of the county. The 
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comparable USAGE Category 5 Immediate Response scenario generates a 26-hour 
clearing time using a network supply of 63,578 vehicles evacuating to outside of the 
county (Mobile District, U.S. Army Corps of Engineers 2010). The MCE-CA clearing 
time for Baldwin County is 77% shorter than USAGE’S clearing time. 

Both MCE-CA estimates are well below the clearing times generated by USAGE. 
However, the MCE-CA’s clearing times represent a utopian scenario in which all 
residents cooperate to minimize the number of stranded residents in each epoch and 
minimize the total travel time of residents. The USAGE model, on the other hand, 
restricts residents to specific evacuation routes. Because MCE-CA has no restrictions on 
resident routing decisions, it generates lower clearing times than USAGE’S restricted 
model. Eor a better comparison, in the next section we explore scenarios in which some 
residents do not evacuate in the most “socially optimal” manner. 

C. IMPACT OF SUB-OPTIMAL ROUTING 

The previous section considered the best-case scenarios for Mobile County, 
Baldwin County, and Escambia County. However, there are several potential reasons 
why best-case clearing times are not achieved during real-life evacuations, such as 
resident response time, resident destination choice and resident routing decisions. In this 
section, we investigate the impact of sub-optimal routing decisions on the part of 
residents. 

I. Scenarios 

Disasters elicit a multitude of evacuation scenarios. According to Cheng et al. 
(2008), residents have many destination choices and may proceed differently in each 
evacuation. Most studies, including the USAGE Mobile County District study of coastal 
Alabama, assign evacuees to destinations using specific routes (Cheng, Wilmot and 
Baker 2008). The USAGE Mobile County District study states that residents tend to 
“follow the leader” during evacuations, unlike normal workday travel when they may 
venture to choose alternate routes (Mobile District, U.S. Army Corps of Engineers 2010). 
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In Section A, we model residents making optimal routing decisions and 
evacuating via any route possible. However, in this section, we model suboptimal routing 
decisions by adjusting MCF-CA’s objective function to incentivize residents to follow 
their own shortest paths rather than the “socially optimal” routing. Let SP^ denote the set 
of arcs (/, j)on the shortest path for resident r, and let w denote the incentive for 
residents to use arcs on their shortest path. Then, the updated MCF-CA objective function 
is: 

Z Z Z (5.1) 

’ ’ (i,J)eA,reR, neN,reR, (iJ)eSPr, 

lsL,teT tsT leLjeT 

Incentivizing residents to use their shortest path effectively decreases the 
flexibility in route choice for residents, similar to studies that use designated evacuation 
routes. We examine clearing times for various incentive levels, with incentive levels 
selected so as to elicit particular routing behavior on the part of the residents. In 
particular, for each county, we compare four scenarios: 

• Global Cooperation: Residents collectively optimize their evacuation, 

• Moderate Cooperation: Approximately 50% of residents use their shortest 
path, 

• Low Cooperation: Approximately 75% of residents use their shortest path, 
and 

• No Cooperation: All residents use their shortest path. 

Note that the Global Cooperation scenario is equivalent to that considered in the 
section A, while the other scenarios consider progressively more myopic routing 
decisions on the part of the residents. 

2. Impact of Sub-optimal Routing on Clearing Times 

Figures 13-15 show that clearing times increase dramatically more residents use 
their shortest path. Escambia County’s clearing time ranges from 7 to 20 hours, Baldwin 
County’s clearing time ranges from 6 to 23 hours, and Mobile County’s clearing time 
ranges from 14 to 48 hours. For all counties, we see that in the No Cooperation scenario, 
clearing time is two to four times that of the Global Cooperation scenario. Low 
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Cooperation decreases clearing time by 6-39% from that of No Cooperation, while 
Moderate Cooperation decreases clearing times by another 13-21%. Note that 
evacuation percentages for all suboptimal routing scenarios quickly diverge from those of 
the Global Cooperation scenario as early as the second hour of evacuation. This finding 
is important because it highlights the important of residents’ route choices in establishing 
clearing times. 


Escambia County 
Evacuation Progression 



•Global Cooperation 
-Moderate Cooperation 
-Low Cooperation 
-No Cooperation 


Figure 13. Escambia County hourly evacuation progression for all scenarios. 


Baldwin County 
Evacuation Progression 



0 5 10 15 20 25 

Hours Elapsed 


—♦—Global Cooperation 
-■-Moderate Cooperation 
-^Low Cooperation 
->^No Cooperation 


Figure 14. Baldwin County hourly evacuation progression for all scenarios. 
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Mobile County 
Evacuation Progression 



0 10 20 30 40 50 

Hours Elapsed 


—♦—Global Cooperation 
-■-Moderate Cooperation 
-^Low Cooperation 
->^No Cooperation 


Figure 15. Mobile County hourly evacuation progression for all scenarios. 


Mobile County’s No Cooperation clearing time grows to more than 200% of its 
Global Cooperation time. Figure 16 reveals that almost 100% of the Mobile County’s 
residents are closest to two of the four designated exit points, I-65N and U.S.-98W. So, as 
more residents follow their shortest path, the remaining two Mobile County exits, U.S.- 
45N and U.S.-43N, are under-utilized. We also observe that Mobile County’s Global 
Cooperation flow use exits more efficiently to start evacuating more residents more 
quickly. By the eighth hour of evacuation. Global Cooperation evacuates 20% more 
residents than Moderate Cooperation. And by hour 14 when Global Cooperation 
completes evacuation, the other scenarios have evacuated just a little over 50% of 
residents. 

This finding is important because it demonstrates that while having many exit 
points established is desirable to having a few, residents must sometimes be encouraged 
to use exits that are distant to ensure timely evacuation of all residents. Thus, intelligent 
routing decisions by residents are the key to ensuring an efficient evacuation. 
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Mobile County Evaucation Exit Use 
All Scenarios 



Global Moderate Low No 

Cooperation Cooperation Cooperation Cooperation 


■ Mobile: I-65N 

■ Mobile: US-98W 

■ Mobile: US-45N 

■ Mobile: US-43N 


Figure 16. Mobile County exit use for all scenarios 


3. Suboptimal Routing Validation and Comparison with Existing 
Models 

We again compare MCF-CA’s clearing time results for Mobile County and 
Baldwin County to those described in the US ACE 2010 study of coastal Alabama 
(Mobile District, U.S. Army Corps of Engineers 2010). We use the USACE 1999 
Northwest Elorida Hurricane Evacuation Study Technical Data Report (U.S. Army Corps 
of Engineers 1999) to validate Escambia County clearing times. We use the clearing 
times of the Global Cooperation and No Cooperation scenarios, shown in Eigure 17, as 
lower and upper bounds for comparison. 
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County Evacuation Clearing Times 


Global Cooperation 

Lower Bound 


No Cooperation 

Upper Bound 




Figure 17. MCF-CA-generated clearing times for all single-county evacuations. 


4. Mobile County Clearing Time 

MCF-CA’s clearing time for Mobile County ranges from 14 to 48 hours. As 
discussed in section B, the comparable US ACE’s Category 5 Immediate Response 
scenario generates a 20-hour clearing time (Mobile District, U.S. Army Corps of 
Engineers 2010), which falls within the upper and lower bounds of the times generated 
by MCE-CA. However, US ACE evacuates only 132,932 vehicles to outside the county, 
while MCE-CA evacuates 197,483. 

We see that US ACE’s Category 3 Immediate Response scenario generates a 
17-hour clearing time and evacuates 10,193 fewer vehicles to outside the county than the 
Category 5 scenario (Mobile District, U.S. Army Corps of Engineers 2010). Thus, we 
assume that the additional clearing time hours are solely attributed to the increased 
number of residents evacuating outside the county. Namely, an increase of 10,193 
vehicles adds 3 hours to Mobile County’s clearing time. Although we cannot extrapolate 
this finding to determine the impact of an additional 64,551 vehicles, we can surmise that 
a complete evacuation of all residents (almost 200,000 vehicles) is well beyond the 20- 
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hour clearing time cited by the USAGE. Moreover, the clearing times generated by MCF- 
CA, 14 to 48 hours, are consistent with the USAGE’S results. 

5. Baldwin Gounty 

MGE-GA’s clearing time for Baldwin Gounty ranges from 6 to 23 hours. 
usage’s comparable Gategory 5 Immediate Response scenario generates a 26-hour 
clearing time (Mobile District, U.S. Army Gorps of Engineers 2010), well above the 
times generated by MGE-GA. Surprisingly, although Mobile Gounty evacuates more than 
twice the number of vehicles as Baldwin Gounty (197,483 and 98,304, respectively), 
usage’s estimated clearing times for Baldwin Gounty are significantly greater than 
those for Mobile Gounty. 

Baldwin Gounty’s USAGE clearing time accounts for normal Mobile Gounty 
through-traffic flow. We test the impact of higher non-evacuation traffic by decreasing 
the capacity of the arcs leading to both of Baldwin Gounty’s exit by 50%. Glearing time 
increases slightly to 14-32 hours when we reduce capacity of arcs leading to I-lOW and 
I-65N. 

While these results are similar to the 26-hour clearing time given by the USAGE, 
additional corroboration is desirable. Unfortunately, the most recent post-hurricane 
assessment for a Gategory 3 or higher storm is the 2009 Post Storm Assessment: 
Hurricanes Gustav and Ike (Mobile District, U.S. Army Gorps of Engineers 2009). The 
report does not provide usable data to assess evacuation clearing times for Baldwin 
Gounty. It states that Baldwin Gounty evacuated “zero” and a “small number” of 
residents for Hurricanes Ike and Gustav, respectively (Mobile District, U.S. Army Gorps 
of Engineers 2009). 

6. Escambia Gounty 

The 1999 USAGE HES for Northwest Florida shows Escambia Gounty’s 
projected 2005 clearing time for a Gategory 5 hurricane is 16.5 hours (U.S. Army Gorps 
of Engineers 1999). The clearing time is based on 100% evacuation, rapid response, and 
low seasonal occupancy. Population growth in addition to road improvements since 1999 
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may cause the clearing times to be slightly different from those reported in the HES. 
However, the most current USAGE HES clearing time is within the 7-20 hour range 
generated by MCE-CA. 

D. OCCURRENCE AND BENEFITS OF CONGESTION 

In this section we examine when and how congestion occurs during the 
evacuation process. Using the four scenarios developed in the previous section, we 
investigate congestion when residents cooperate and when they do not. We then restrict 
MCE-CA to disallow congestion in all arcs in all time epochs. Restricting MCE-CA to 
disallow congestion models an ideal controlled departure system that maximizes vehicle 
travel speed during evacuations. We study the impact of this restriction on clearing times. 

I. Impact of Congestion with Sub-optimal Routing 

Eigure 18 displays Mobile County residents’ evacuation progression and the 
percentage of roads utilized in each epoch that are road congested. We see that 
congestion persists throughout the evacuation horizon. We expect to see higher levels of 
congestion associated with long clearing times and low levels of congestion associated 
with short clearing times. However, that assumption is incorrect for all county 
evacuations. 

The percentage of utilized roads that are congested for both the Low Cooperation 
and Moderate Cooperation scenarios rise above No Cooperation at various points in the 
evacuation process. The first time Low Cooperation rises above No Cooperation is in 
hour 12. Up to this point, approximately 47% of residents have evacuated in both 
scenarios. At hour 12, additional residents start moving on the road network in Low 
Cooperation while no additional residents start moving in the No Cooperation. More 
roads are congested in the Low Cooperation scenario simply because more arcs are 
utilized. We see similar results for both Escambia County and Baldwin County (see the 
Appendix). 
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Mobile County 

Global Cooperation Moderate Cooperation ^^^Low Cooperation ^^No Cooperation 


Evacuation Progression 



Road Congestion 



Figure 18. Mobile County evacuation progression and road congestions. 

2. Resident Congestion 

While it is useful to exa mi ne the percentage of road segments that experience 
congestion, it is perhaps more instructive to consider the percentage of residents that 
encounter congestion during their evacuation. Figures 19-21 show the percentage of 
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residents in each county that experience congestion at some point during the evacuation 
horizon and the total number of arcs that are congested. We see that more than 50% of 
residents in each county experience congestion. 

We observe that the percentage of residents who experience congestion is lowest 
in scenarios in which more arcs are congested. Global Cooperation has the lowest 
percentage of residents that experience congestion, the shortest clearing time, and the 
highest number of congested arcs. The flexibility in route choice that Global Cooperation 
affords increases the number of arcs used, thereby increasing the number of congestible 
arcs. This result implies that increased flexibility in route choice improves clearing time. 


Baldwin County 


Residents Experiencing Congestion 



Cooperation Cooperation Cooperation Cooperation 
(221 ares) (107 ares) (89 arcs) (78 arcs) 


Figure 19. Baldwin County’s proportion of residents that experience congestion 

and the number of congested arcs. 
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Escambia County 


Residents Experiencing Congestion 



Cooperation Cooperation Cooperation Cooperation 
(319 ares) (223 ares) (134 arcs) (124 arcs) 


Figure 20. Escambia County proportion of residents that experience congestion 

and the number of congested arcs. 



Figure 21. Mobile County proportion of residents that experience congestion 

and the number of congested arcs. 
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3. Benefits of Congestion on Clearing Times 

A common focus of evacuation research is on dispensing evacuation notices in 
such a way as to minimize congestion (Li et al. 2010). In this section, we exa mi ne 
whether clearing time is worsened by congestion, or improved. We accomplish this by 
restricting MCF-CA to disallow congestion on all arcs in all time epochs, i.e., by setting 
=1 for all We denote the resulting clearing times as the non-congested 

clearing times and compare them to the congested clearing times generated in the 
previous section. 

Figure 25 depicts congested and non-congested clearing times during a mass 
evacuation for all scenarios. Clearing times for the non-congested evacuations range 
between 22 and 130 hours, while clearing times for the congested evacuations range 
between 15 and 83 hours. By the fifth hour of evacuation, congestion helps evacuate 
22-43% of residents, while only 16-28% of residents evacuate in the non-congested 
scenarios. 

One might postulate that the 70-80% decrease in travel speed when going from 
under- to over-saturated conditions would sharply increase the congested clearing time. 
However, in over-saturated conditions over 50% more residents are able to start or 
continue evacuating in each time period, albeit at a slower speed. This finding is 
important because it reveals that while it is an annoyance, congestion alone does not 
prolong clearing times beyond those of non-congested evacuations. In fact, congestion 
improves evacuation clearing times by 20-36%. Therefore, even if state and local 
emergency managers were able to design a communication system that could produce 
100% resident compliance and alleviate congestion, they would need to carefully weigh 
the impact on clearing time before implementing such a system. 
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Mass Evacuation Progression 
Congestion vs. No Congestion 



Hours Elapsed 


NO CONGESTION 

- Global Cooperation 

- Moderate Cooperation 

- Low Cooperation 

- No Cooperation 
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— • -Global Cooperation 
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— • -Low Cooperation 

— • -No Cooperation 


Figure 22. Evacuation Progression for congested and non-congested mass evacuations. 
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VI. CONCLUSION AND FUTURE WORK 


State and local emergency managers study evacuation clearing times in order to 
establish appropriate traffic control measures to mitigate the congestion that is likely to 
occur during the mass movement of residents. This thesis develops a new model for 
evacuations and uses this model to investigate evacuation clearing times on the multi- 
regional transportation network of Mobile County, AL, Baldwin County, AL, and 
Escambia County, FL. We create a spatial representation of this transportation network 
and compute the number of vehicles using this network using household vehicle 
availability data from the 2010 Census Survey. 

We develop a multi-period integer linear program called the Minimum Cost Flow 
with Congestion Scheduling (MCF-CA) model. MCF-CA discretizes the nonlinear 
relationship between traffic flow, speed, and density. MCF-CA selects the lowest optimal 
congestion level that will maximize flows (i.e., throughput) in order to minimize the 
number of residents not evacuated for each time period. 

We measure the impact of mass evacuation on county clearing times and explore 
the impact of congestion and routing on clearing times. Finally, we validate results using 
the 2010 clearing times generated by the U.S. Army Corps of Engineers, Mobile County 
District. 

A. SUMMARY AND RECOMMENDATIONS 

During evacuations, mass or otherwise, congestion is inevitable. In the best- and 
worst-case scenarios, more than 50% of residents will be congested at some point during 
the evacuation horizon. Nevertheless, this thesis demonstrates that congestion alone does 
not increase clearing times. Clearing times are shorter when residents are allowed to 
evacuate earlier when congestion is allowed versus later when congestion is not allowed. 
Congesting improves evacuation clearing time by 20-36% over not congesting because 
50% more residents are able to start or continue evacuating during each time epoch. 
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However, MCF-CA does not consider secondary effects such as car accidents due 
to closer vehicle spacing during congestion. States and counties have many options to 
handle congestion, including improving roads, constructing new corridors, and using 
contra-flow. However, due to the growing rate of congestion in cities, new facilities 
would soon face the same congestion. Therefore, to keep the time-saving benefit of 
congesting, emergency managers might study traffic control measures that would help 
mitigate the secondary effects of congestion. 

We observe that exit locations, not simply the number of exits, play a significant 
factor in estimating clearing times. Escambia County’s clearing time during mass 
evacuations rise as much as threefold because a third of its exit points border Baldwin 
County. Having one set of exit points for both county and mass evacuations would 
decrease the difference in clearing times. Additionally, standardizing for both evacuation 
types would make it easier for residents to remember evacuation routing and procedures 
for any scenario. 

Finally, we show that in scenarios with a higher percentage of residents following 
their directed shortest path routes, clearing times lengthen by as much as two to four 
times. Residents must sometimes be encouraged to use exits that are distant in order to 
ensure timely evacuation of all residents. Thus, intelligent routing decisions by residents 
are the key to ensuring an efficient evacuation. While it is not realistic to expect that 
residents will cooperate as a system to minimize clearing times, it is equally unrealistic to 
assume that residents will stay on a designated route throughout the evacuation horizon. 
As a compromise, MCF-CA could be implemented in future studies to analyze the effect 
of varying the amount of residents who decide not to follow the designated route to 
specified destination points. 
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B. FUTURE WORK 


The MCF-CA model has enough flexibility to be applied on all levels of analysis 
(microscopic, mesoscopic, and macroscopic). Some areas of possible research are: 

• Determining the impact of individual route-change decisions when faced 
with traffic congestion, 

• Exploring the effect of individual household destination decisions on 
traffic congestion, 

• Analyzing the effect of varying the number of individuals who follow 
designated routes to specified destinations, and 

• Designing a community controlled departure system that decides route and 
departure time jointly and accounts for disaster uncertainty. 
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APPENDIX. ADDITIONAL GRAPHS 


C. CONGESTION DURING SUB-OPTIMAL ROUTING 

Figures 23 and 24 show Baldwin County and Escambia County residents’ 
evacuation progression and the percent of roads utilized in each epoch that are road 
congested, respectively. 


Baldwin County 

■Global Cooperation M Moderate Cooperation A Low Cooperation )( No Cooperation 


Evacuation Progression 



Road Congestion 



Figure 23. Baldwin County evacuation progression and road congestions. 
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Escambia County 

■Global Cooperation M Moderate Cooperation A Low Cooperation )( No Cooperation 


Evacuation Progression 



Road Congestion 



Figure 24. Escambia County evacuation progression and road congestions. 
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